Three-body rf association of Efimov trimers 
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We present a theoretical analysis of rf association of Efimov trimers in a 2-component Bose gas with 
short-range interactions. Using the adiabatic hyperspherical Green's function formalism to solve the quan- 
tum 3-body problem, we obtain universal expressions for 3-body rf association rates as a function of the 
s-wave scattering length a. We find that the association rates scale as a~ 2 in the limit of large a, and 
diverge as a ajj^ whenever an Efimov state crosses the atom-dimer threshold (where a ad stands for the 
atom-dimer scattering length). Our calculations show that trimer formation rates as large as ~ 10 -21 cm 6 /s 
can be achieved with rf Rabi frequencies of order 1 MHz, suggesting that direct rf association is a powerful 
tool of making and probing few-body quantum states in ultracold atomic gases. 
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The experimental study of universal phenomena in ul- 
tracold atomic gases [Q], y] has led to numerous recent ad- 
vances in quantum few -body physics llfl. In particular, 
weakly bound Efimov trimers formed by resonantly inter- 
acting bosons J2] or fermions yt] provide a unified frame- 
work for understanding 3-body phenomena in atomic, nu- 
clear, condensed-matter, and high-energy physics. The 
field of ultracold quantum gases provides an ideal play- 
ground for studying the Efimov physics because interac- 
tions of ultracold atoms can be tuned over many orders 
of magnitude by varying an external magnetic field near 
a Feshbach resonance H. Measurements of 3-body re- 
combination rates as functions of the atom-atom scatter- 
ing length a have provided solid evidence for the existence 
and universal properties of the Efimov states in ultracold 
gases of alkali-metal atoms J2, @-@], providing novel in- 
sights into the quantum dynamics of 3- and 4-body recom- 
bination 0]. 

While highly successful, experimental studies of Efimov 
trimers via 3-body loss measurements are limited to ob- 
servations of zero-energy crossings of the trimers with the 
continuum J2,@, El- Very recently, this limitation has been 
overcome by Lompe et al. and later by Nakajima et al. 
JH], who reported signatures of Efimov states in radiofre- 
quency (rf) association spectra of an ultracold mixture of 
6 Li 2 dimers with 6 Li atoms. These experiments have pi- 
oneered a radically new approach to few-body physics, in 
which Efimov states are "actively" created and probed via 
rf spectroscopy flU, rather than "passively" observed via 
measurements of 3-body recombination rates. This ap- 
proach is potentially very powerful, since rf spectroscopy 
is a well-established experimental tool for probing few- 
and many-body interactions in ultracold quantum gases 
U [l(|. Examples include measurements of binding en- 
ergies of Feshbach dimers using bound-bound and bound- 
free rf spectroscopy fl, as well as many -body correlations 
and quantum phase transitions in ultracold quantum gases 
iflQfl - Near-resonant rf fields can be used to create Fes- 



hbach molecules by rf association 111 111 , and manipulate 
magnetic Feshbach resonances in collisions of ultracold 
atoms il2lll3ll . 

In this Letter, we develop a quantum theory of 3-body 
rf association of Efimov trimers. By combining the adia- 
batic hyperspherical Green's function formalism [fl4ll with 
Fermi's Golden Rule, we derive the expressions for 3-body 
association rates as a function of the atom-atom scattering 
length a. Our results show that 3-body rf association rate 
scales as a~ 2 in the limit of large a and approaches infinity 
whenever an Efimov state crosses the atom-dimer thresh- 
old at positive a. We apply our theory to calculate the rate 
for the production of Efimov trimers by rf association in 
an ultracold gas of 7 Li atoms. Our results show that large 
rf association rates can be achieved with Rabi frequencies 
available in the laboratory. Our theory can be generalized 
to describe rf-assisted 3-body phenomena in 3-component 
Fermi gases and atom-molecule mixtures [SO], providing 
a novel tool for studying rf association, spectroscopy, and 
chemical dynamics of universal quantum few-body states 
in ultracold atomic gases. 

We begin by considering rf association of three free 
atoms b+b + x leading to the formation of a weakly bound 
Efimov trimer bbb where b and x refer to two different hy- 
perfine states of a bosonic atom. The experimental proto- 
type for such a system is a 2-component Bose gas of 7 Li 
atoms in their lowest hyperfine states 0, [H]. The system 
of three particles is characterized by two s-wave scattering 
lengths abb an d o-bx- Throughout this work, we assume that 
o-bx = 0, and examine the dependence of the scattering 
observables on a single parameter - the scattering length 
a = a b b- The thermal rate constant for rf association of 
three distinguishable atoms into a weakly bound trimer can 
be written in the form 

hk hk 32-7T 2 . _ . „ A . ., 
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where o { f is the cross section for stimulated radiative as- 
sociation, [i is the 3-body reduced mass, k 2 = 2fj,e/h 2 , e 
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is the collision energy, 0, the Rabi frequency of the rf field, 
A = Hto — €f is the detuning from resonance, and ef is 
the binding energy of the Efimov trimer formed by rf as- 
sociation. While £Q is similar in form to the expressions 
encountered in the theory of photoassociation lineshapes 
1 17, 18h . the important pref actor 32n 2 /k 5 have been intro- 
duced to properly account for the 3-body phase space den- 
sity lfl9h . In deriving Eq. £Ql, we assumed that the temper- 
ature of the ultracold gas is sufficiently low so that only s- 
wave collisions contribute to the cross section. It is conve- 
nient to parametrize the S'-matrix element by a Lorentzian 
[TSUI 
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where 7«(e, f2) is the stimulated width, 7< j is the decay 
width introduced to account for non-radiative decay of 
Efimov trimers via 3-body recombination [21], and 7 = 
7s + Id is the total width. The stimulated width is given by 
the Fermi's Golden Rule (vA\l$l 

ls (e,n) = 2n^-\^ l \^ f )\ 2 (3) 



where = 2{x\jl\b) /h is the atomic Rabi frequency, /t 
is the transition dipole moment operator, ^ \ is the wave- 
function for the initial 3-body continuum state, and *S> f 
is the wavefunction of the Efimov trimer. In the limit of 
zero collision energy, | Sif\ 2 ~ e 2 ||20|| . and we can define 
an energy-independent normalized Franck-Condon factor 
(FCF) 



T if = \(V i \*f)\ 2 /e 2 



(4) 



Substituting this expression in Eq. CO, we obtain for the 
3-body association rate on resonance (e = A) 



(5) 



M° Id 
In deriving this expression, we made the assumption 7^ <C 
7d, which is well justified for short-lived Efimov states 
with observed decay widths on the order of several MHz 
flm. In order to evaluate the FCF ©, we use the adia- 
batic hyperspherical Green's function formalism IU4II . We 
assume that the initial and final 3-body wavefunctions take 
the form * p = R~ 5 / 2 F p (R)^ p (R]to), where R is the 
hyperradius, to are the hyperangles, and (3 = i,f. Here 
$i is assumed to the hyperangular channel function which 
corresponds to three free particles in the R^§> a limit, and 
$ f is taken to be the channel function which supports an 
Efimov trimer state. Substituting these wavefunctions to 
Eq. ©, we find 
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where Wi f (R) is the hyperangular channel function over- 
lap at a fixed hyperradius and the radial functions Fa (R) 
are solutions to the adiabatic Schrodinger equations Qj]] 

U f3 (R))F fj (R) = e p F f3 (R), (7) 



2fi \dR 



subject to the boundary condition Fp(Ro) = 0, where R 
is the radial cutoff parameter related to the experimentally 
relevant 3-body parameter [HI by R ~ 0.22a + where 
a + is the position of the first universal minimum in the 
3-body recombination rate of three identical bosons. At 
R — > 00, the energy-normalized scattering wavefunction 
takes the form F t (R) ~ (2 / u/vrft 2 A:) 1 / 2 sm(kR + 5) §], 
where 5 is the scattering phase shift. For the bound-state 
wavefunction, we impose the standard boundary condition 
^(00) = 0. In Eq. @, Up{R) = [(v p (R) + 2) 2 - 
1/4] /R 2 + Q(R) are the adiabatic hyperspherical poten- 
tials, vp{R) are the adiabatic eigenvalues and Q{R) are 
the diagonal non-adiabatic corrections II 14-11 . The angular 
integral Wif in Eq. (O can be evaluated by expanding the 
functions ^^(R; to) in hyperspherical harmonics llll[l4ll. 

Figure [T] shows the calculated FCFs for the first two 
Efimov states in unites of the dimer binding energy e^ 3 
(e D = h 2 /ma 2 ) as a function of the atom-atom scattering 
length in units of the 3-body cutoff R a . As a/R increases, 
the Efimov states become more strongly bound, and the 
bound-state wavefunction ^ f(R) in Eq. (O becomes more 
localized in the small R region, where the amplitude of the 
scattering wavefunction *$>i(R) is suppressed due to a steep 
increase of the adiabatic potential Ui(R) lfl4ll . The angu- 
lar overlap integral Wif(R) becomes constant in the limit 
a — > 00 (see the inset in Fig. (TJ, so the large-a behavior of 
the FC overlap is determined solely by the wavefunctions 
in Eq. (©. The decrease of the FCF with a shown in Fig. Q] 
is therefore a consequence of diminishing overlap between 
the initial and final wavefunctions. 

The FCFs for different Efimov states have similar 
shapes, decreasing monotonously with increasing a/R Q . 
The FCF for the first excited Efimov state is shifted to- 
ward higher a by the universal factor of 22.7 [!!]. Since 
the FCFs scale as e^ 3 we can expect that the FCFs for dif- 
ferent Efimov states will be related by a factor of 22. 7 6 . 
However, in general the decay width of consecutive Efimov 
states will decrease by the universal scale factor of 22. 7 2 
meaning that the zero temperature 3-body association rate 
constant in Eq. (|5]l will scale as 22. 7 4 between neighboring 
Efimov states. 

To make these arguments more quantitative, let us first 
consider the limit of small a, where the trimer state ap- 
proaches the atom-dimer threshold and the radial wave- 
function is well approximated by 



F f (R) ~ \J 4^ (e, - e D ) exp 2/i (e f - e D )R 

(8) 

In this limit, most of the overlap occurs when R S> a and 
thus we approximate the initial energy-normalized wave- 

nr 

function as F{ (R) 



RJ2 (kR), where J2 is a 



Bessel function. Examining the numerical results shown 
in the inset of Fig. ffl yields W (R) -> Ca 3 / 2 /R 3 ^ 2 with 
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Fig. 1. The FCF in units of e^, 3 as a function of the reduced 
scattering length a/Ro- Symbols - numerical calculations using 
Eq. ©; dashed lines - analytic results obtained for a — »■ (blue 
dashed line) and a — > oo (magenta dot-dashed line). The inset 
shows the radial dependence of the angular integral Wif(R) in 
reduced units. 
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452.741. Using this result in Eq. ©, we find 

2e + e-2Ve (e + e))' 
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where e = 6f — is the binding energy of the Efimov 
state relative to the atom-dimer threshold. The analytical 
result shown in Fig. [T] is in excellent agreement with the 
full numerical calculation; the difference does not exceed 
10% for a/R < 2.5. In the limit where ka <C 1 and 
e <C e Eq. (O takes on the simple universal form Ti$ — > 

^ 3 C 2 a 3 a 3 ad / (V23 3 / 4 ft 6 ) where a ad = y/3H 2 /Am£ is 

the scattering length of a bosonic atom scattering of a 
boson-boson Feshbach molecule. As the first Efimov state 
crosses the atom-dimer threshold and becomes unbound at 
a/R ~ 2.5, a ad tends to infinity and J 7 ^ diverges as 
shown in Fig. [TJ 

We now consider the opposite limit a — > oo, in which 
the Efimov state becomes strongly bound relative to the 
atom-dimer threshold and the final wavefunction is only 
non-zero in the region R < a. Further, the adiabatic po- 
tentials in Eq. become universal when R <C a, giving 
the initial and final wavefunctions 

n 2 




Fi(R) 



(10) 

where K iso is a modified Bessel function of the second 
kind of imaginary order and s 2 = 1.01251 flUl- As 



shown in the inset of Fig. [TJ W (R) ->• B v with B x = 
47.898 in the limit a — > oo. Substitution of this result in 
Eq. © gives 

Bl ( eV sinh (vrso) T 2 T 2 (u + ) 
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where v± = (v + 4 ± is )/2 and 2 Fi is a hypergeometric 
function. This expression for the FC overlap is correct in 
the case where \ka\ 3> 1 not only for the v = 1 initial 
channel of interest in this work, but also for all incident s- 
wave channels with v = —1, 1, 3, 5, .... By replacing €f 

with €f n = £/ (22. 7) 2 ™ we can get the overlap of the nth 
Efimov state in terms of the first with binding ef. In the 
large a limit, but where ka <C 1, there is extra suppres- 
sion of the wavefunction in the u = 1 incident channel due 
to the additional tunneling in the a < R < 1/k regime 
producing and extra factor of (6/ka) 2 in the FCF. When 
a — > oo and in the ka <C 1 regime the FC overlap takes 
on the simple universal form of J 7 ^ — > 1661. 127l 2 / \xa 2 € A ^. 
It should be noted that this large a behavior was found un- 
der the assumption that the initial 3-body state is described 
by a single hyperangular channel function. In truth, in the 
region where R <C a the initial state of the system is multi- 
channel in nature, and a more complete description will be 
required. Fortunately the FCFs found in Eq. ([Til are fully 
generalized for overlaps in this regime. Incorporating the 
full multi-channel nature of the system is the subject of on- 
going study. 

As an application of the universal theory developed 
above, we evaluate the rate for the production of Efimov 
trimers by rf association in an ultracold gas of 7 Li atoms. 
The Efimov physics in this system has been studied in de- 
tail experimentally flBh. We assume that the 7 Li atoms 
are prepared in two lowest hyperfine states \b) = |11) and 
\x) = 1 10) and the dc magnetic field is tuned such that 
o-bx = 0. Figure [2a) shows the 3-body rf association rate 
as a function of a calculated using typical experimental pa- 
rameters Q/27T = 10 kHz y|], Ro = 58ao, and 7^ = 4.5 
MHz fl flal . The association rate takes a minimum value 
of 3.3 x 10~ 26 cm 3 /s at a = 200a and a maximum value 
of 2.9 x 10~ 25 cm 3 /s at a = 1940a - These rate constants 
are large enough as to be easily observable experimentally 
via monitoring atom loss from an optical dipole trap Q4|] . 
The oscillating behavior of K 3 observed at finite a is in 
marked contrast with the familiar o 4 scaling of the 3-body 
recombination rate |l]|. 

As the trimer binding energy approaches the atom-dimer 
threshold, it is natural to expect that rf association will start 
producing not only Efimov trimers, but also weakly bound 
dimers. In order to examine the interplay between these 
processes, we evaluated the 2-body rf association rate K 2 
using the Fermi's Golden Rule iHdUl in the limit b<l 
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Fig. 2. (a) The 3-body if association rate for Li as a function of 
a calculated for Q/2ir = 10 kHz and R = 58a H|. (b) The 
relative efficiency of 3-body vs 2-body if association (K-znjKi) 
plotted as a function of a for n = 10 12 cm~ 3 and T = 1/iK. 



1 15] with the result 



K 2 



(12) 
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In Fig. [2b), we plot the ratio g = K 3 n/K 2 , which quan- 
tifies the relative efficiency of 2-body versus 3-body rf as- 
sociation, for typical experimental conditions n = 10 12 
cm~ 3 and T = 1 fiK |4J] (note that g is independent of 
fi). The ratio declines monotonously with a due to the a 3 
scaling behavior of the 2-body association rate (fT2t . The 
results shown in Fig. |2]thus suggest that 3-body rf associa- 
tion experiments should be performed at small-to-moderate 
values of a to suppress 2-body rf association. 

In summary, we have developed a universal theory of rf 
association of Efimov trimers in a 2-component Bose gas 
with short-range interactions. We have applied the theory 
to evaluate the rates for 3-body rf association in an ultra- 
cold gas of 7 Li atoms and derived analytical expressions 
for the association rates, which are valid in the limits of 
large a (K ~ a -2 ) or infinitely weakly bound Efimov 
states (K ~ a 3 a^ d ). We find that rf association of weakly 
bound (atom-dimer-like) Efimov trimers is more efficient 
than that of strongly bound Efimov trimers (Fig. |2j, and 
rationalize this result based on simple overlap arguments. 
We also examine the competing process of 2-body rf as- 
sociation, and find that the relative efficiency of 3-body 
versus 2-body rf association decreases with increasing a 
(Fig. [2J})- The methods developed in this work can easily 
be generalized to describe rf-induced few-body phenomena 



in 3-component Fermigases JH], atom-dimer mixtures 1161 
and 4-body systems |7|] or in 3-body systems with multi- 
channel 2-body interactions Jill!. These methods might 
also be applied to isotropically trapped systems where a 
transition matrix element to an excited trap state can be 
calculated in stead of a rate constant. 



Our findings provide quantitative insight into the mech- 
anisms of rf association of universal trimer states in ul- 
tracold quantum gases, and may thus have significant im- 
plications for research in this rapidly expanding area of 
physics IB-Ill]. We have shown that rf association of ul- 
tracold Li atoms into Efimov trimers occurs at a substan- 
tial rate of ~10 -25 cm 6 /s (Fig. |2]l for Rabi frequencies of 
order 10 kHz, which can be readily achieved in the labo- 
ratory JUdHl. The 3-body association rate may be further 
increased by raising the Rabi frequency, since K 3 ~ Q 2 
(0, so values of K 3 on the order ~ 10~ 21 cm 6 /s should be 
readily achievable with Q/2ir ~ 1 MHz. Such Rabi fre- 



quencies can be generated in atom chip experiments 111211 . 
which would provide an ideal setting for creating large 
quantities of Efimov trimers, heretofore an unsolved prob- 
lem due to their large decay rates yfl. The results shown in 
Fig. |2] suggest that rf association experiments in Li should 
be performed at small-to-moderate values of a (between 
200 and 2000a ) to avoid the competing processes of 2- 
body rf association and 3-body recombination. 
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Smithsonian Astrophysical Observatory. 
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